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We analyze the singlet-triplet splitting on a planar array of quantum dots coupled capacitively 
to a set of external voltage gates. The system is modelled using an extended Hubbard Hamiltonian 
keeping two excess electrons on the array. The voltage dependence of the low-energy singlet and 
triplet states is analyzed using the Feshbach formalism. The formation of a well decoupled two- 
level system in the ground state is shown to rely on the fact of having two particles in the system. 
Coherent operation of the array is studied with respect to single quantum bit operations. One 
quantum gate is implemented via voltage controls, while for the necessary second quantum gate, 
a uniform external magnetic field is introduced. The Aharonov-Bohm phases on the closed loop 
tunnel connections in the array are used to effectively suppress the tunneling, despite a constant 
tunneling amplitude in the structure. This allows one to completely stall the qubit in any arbitrary 
quantum superposition, providing full control of this interesting quantum system. 



PACS numbers: 03.67.Lx, 03.65. Vf, 85. 35. Be, 85.35.Gv 



I. INTRODUCTION 

Semiconductor quantum dot systems have been stud- 
ied extensively in recent years. These quantum dots 
(qudots) typically contain from a few to a few hundreds 
of extra-electrons and are externally controlled by volt- 
age gates. Their behavior at low bias voltages can be 
well understood by looking at the energetically topmost 
electrons in the qudot i For an even number of electrons, 
they may pair up such that the total spin is zero and the 
topmost two electrons essentially form a singlet, a system 
that has been proposed as a possible source for entan- 
gled electrons. 2 However, under special circumstances, 
the ground state may actually be a triplet even without 
the presence of an external magnetic fields A straightfor- 
ward explanation for this behavior can be given in terms 
of the exchange contribution to the energy from the two 
topmost electrons for the case of nearly degenerate single 
particle levelsi For an odd total number of electrons on 
a single qudot, the overall spin is typically found to be 
1/2. 

In this paper we consider an ensemble of qudots, 
specifically a 2 x 2 array of interconnected and interact- 
ing qudots similar to the typical cellular automata unit 
cell format The interaction is modelled within a 

capacitance matrix formalism. In the weak tunneling 
regime, the direct exchange energy from the two-body in- 
teraction between electrons on different qudots is negligi- 
ble. Residence of these two electrons on the same qudot 
(double occupancy) is energetically unfavorable due to 
the comparatively large Coulomb charging energies. The 
2x2 array is considered as a particular implementation of 
a charge quantum bit (charge qubit) . As shown in [t| , it is 
essential for the single qubit operations of such a charge 
qubit to have a non-local potential, tuneable tunneling 
amplitudes or an external magnetic field which provides 
a complex phase to the wavefunction. Since the direct 
exchange contribution with charges on different dots is 
negligible, there are no non-local potential effects here. 



Further, the tunneling amplitude is considered constant, 
fixed by the specific realization of the solid state qudot 
array which typically uses oxide barriers between qudots, 
and are then basically unaltered by potential gates. In 
this context of a fixed geometry, the only way to im- 
plement full single qubit operations is using an external 
magnetic field which can be controlled at willi& It is fur- 
ther considered that the field is uniform within the area 
of the qudot array as is likely the case in normal imple- 
mentations. 

Due to the geometrical symmetry of the qudot array, 
the ground state turns out to be exactly degenerate for 
the triplet states when no magnetic field and no gate 
voltages are applied. However, the spectrum exhibits a 
gap for the singlet states which is related to the distinct 
symmetry under particle exchange for the singlet and 
triplet states. Therefore, particle exchange does play an 
essential role, yet it is a higher order effect, similar to 
superexchangeiS, as the exchange of two particles on the 
qudot array takes more than one tunneling step and ex- 
plores virtual higher energy states. 

The properties of the singlet and triplet states are de- 
termined by two remaining parameters, the asymmetri- 
cally applied gate voltage V g , which breaks the 90 degree 
symmetry, and the external magnetic field B perpendic- 
ular to the array. We show that careful control of the 
magnetic flux through the system allows one to rotate 
the qubit Bloch vector about an axis orthogonal to the 
one provided by V g . These two parameters then, provide 
full quantum manipulation of the charge qubits. 

The analysis of our numerical results follows the Fes- 
hbach formalism (App. ^ which provides an effective 
Hamiltonian of the two relevant dimensions of the qubit 
subspace, well decoupled energetically from higher lying 
states. This approach provides insights on the nature of 
the system dependence on fields, and explains the abil- 
ity of the magnetic field to complete the set of necessary 
single qubit operations. 
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The Hamiltonian used to describe this system is of the 
extended Hubbard type 




FIG. 1: Setup of 2 x 2 array, (a) Arrangement of the four is- 
lands with mutual tunnel connection indicated by black lines, 
(b) Same as in (a) but shows explicitly tunneling junctions 
and capacitive coupling including the two voltage gates act- 
ing along the square diagonals. 



II. THE MODEL SYSTEM 

We examined a range of different geometries of pla- 
nar arrays with capacitively coupled quantum dots. The 
best decoupled quantum 2 level systems (qu2LS) in the 
ground state were found to exhibit high spatial symme- 
try, i.e. Ci V symmetry, in agreement with the require- 
ment of spatially distinct wave functions needed to par- 
ticipate in the qu2LSi^ The Ca v symmetry of the 2x2 
array, for example, ensures a binary groundstate system 
that for weak inter-dot tunneling is well decoupled en- 
ergetically from the remainder of the spectrum. Other 
geometries are clearly possible, but we have chosen here 
one of the simplest closed loops with C± v symmetry to 
demonstrate the main concepts. 

The model network under consideration is a 2 x 2 array 
of qudots with a single spatial state per site plus spin. 
The system is sketched in Fig. QJi. As explicitly out- 
lined in panel (b) , tunneling is allowed between any pair 
of dots, where every tunnel junction also carries capac- 
itance. The parameters which enter the model are: the 
dot-gate capacitance (C g — 25aF), the dot-dot capaci- 
tances (Cdd — 25aF for nearest neighbor dots and 17aF 
for dots connected through the diagonal of the array), the 
dot self-capacitance (Cdo = 25aF) and the dot-dot tun- 
neling amplitude (t = 2fieV) . The parameters have been 
chosen such that the energy cost for double occupancy of 
a dot (standard Hubbard U) is about lmeV, a commonly 
used value for these systems, and corresponding to typ- 
ical dot dimensions of lOOnm. Further, the tunneling t 
must not be chosen too large since otherwise the quality 
of the 2-level system in the groundstate as well as the 
well-defined charge per dot are compromised, as we will 
see below. 



i.{T idt& 

2 H VijC+cfjCj^Cia+^TVifii (1) 

with hi = cXcii + cjcjj and cf a the typical creation op- 
erator for a particle at qudot i with spin a. £u\ a refers 
to the local energy of state a on the i = {1, . . . , n} iden- 
tical dots and can be used to account for the Zeeman 
splitting of spins in an external magnetic field. Through- 
out this paper, however, the e a are simply set equal and 
zero. The tunneling coefficients t^- from dot i to dot j 
are considered independent of the spin orientation, thus 
ij- = t\- = tij. Furthermore, the tunnel connections are 
considered to be the same up to a phase, i.e. |ty| = \t\. 
The electrostatic energy in the last two terms of Eq. QJ, 
i.e. the coefficients Vij and Vi, are derived from the total 
capacitance matrix of the system which is approximated 
by the capacitor network indicated in Fig. \I]p& 

An essential property of singlet and triplet states is the 
effective separation of the spin degree of freedom from 
the spatial wavefunction component since the total state 
is a product of spatial and spin components. Here the 
(anti)symmetry of the spatial wavefunction of (triplet) 
singlet states under particle exchange is taken care of 
by imposing (anti)commutator relations on the creation 
and annihilation operators cf and Ci with the spin in- 
dex dropped, yet keeping the constraint of a total of two 
electrons. 

An external magnetic field perpendicular to the net- 
work of qudots affects the system insofar as spatial prop- 
agation is associated with the acquisition of a complex 
phase. Therefore the become comple»ii 

Uj = t*i = \tij I j*** , with Vij = ^J A -di, (2) 

Xi 

where A is the vector potential of the applied magnetic 
field, B = V x A. Using a symmetric gauge, the acquisi- 
tion of phase in the 2x2 array of qudots is indicated in 
Fig. The phase flows clockwise on the outer connec- 
tions while the diagonal connections remain phaseless. 
The ring structure leads to an Aharonov-Bohm phase 
(AB) for a single particle moving around the ring. Note, 
however, that the second particle is essential for the nec- 
essary ground state qu2LS needed in the qubit setup. For 
lithographic setups on the scale of 200nm, the required 
magnetic field for an AB phase cycle on the whole ar- 
ray is around lOOmT and thus rather small. At these 
fields the local wave functions in the individual quantum 
dots do not change much, and the absolute value of the 
tunneling \t\ is considered constant. 
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FIG. 2: 2 x 2 qudot array, (a) Array with a perpendicular 
magnetic field applied, (b) Schematic network of the Hilbert 
space states with tunneling transitions between them indi- 
cated by connecting lines (states of double occupancy are not 
included). The arrows indicate the flow of complex phase 
acquired from a magnetic field in the tunnelling coefficients 
\t\e lip . The dashed blue lines indicate paths of particle ex- 
change (see text). |0) 6 = 1 13) and \l)„ b = 1 24) indicate the 
qubit states where \ij) are the two electron states with one 
electron on dot i and the other on dot j. 

In Fig. 03, the two-particle Hilbert space and the al- 
lowed tunnel transitions are shown for the 2x2 system, 
for simplicity ignoring states of double occupancy. The 
blue dashed lines in Fig. ^ are related to particle ex- 
change in the sense that the off-diagonal element in the 
corresponding triplet Hamiltonian has an extra minus 
sign due to its fermionic character, as it can also be seen 
directly from the basis chosen. For example, the transi- 
tion 1 12) to 1 23) can be thought of as the two-step process 
of one hopping and one exchange, |12) — > 1 32) — » — |23). 
Any path with an odd number of dashed segments in the 
Hilbert space of Fig. ^> has an extra minus sign associ- 
ated with it. Notice also, that if the path in the Hilbert 
space network of Fig. |2h is closed, then an odd number 
of dashed segments refers to an effective exchange of the 
two electrons. More on this later. 



III. ANALYSIS 

Using the Hamiltonian in Eq. QJ, the eigenspectrum 
for the 2x2 array is shown in Fig. as function of an 
asymmetric gate voltage drag V g , with no magnetic field 
applied (B = 0). For small V g , the singlet (triplet) low- 
energy level set of interest is well separated from the re- 
maining singlet (triplet) spectrum. In contrast to the sin- 
glet set, however, which has an anticrossing at V g = 0, the 
triplet set is degenerate there. Panel (b) shows the eigen- 
spectrum as function of the magnetic field when there 
are no gate voltages applied (V g — 0). The tunneling 
has been chosen relatively large (\t\ — 5/xeV) such that 
the energy splitting S due to the tunneling in the low- 
energy singlet set (the qu2LS we will focus on) reaches 
about 1/10 of the distance to the nearest higher lying 
states, A. This is still a good qubit configuration, as the 
coherent state manipulation in the qu2LS can be per- 



formed without significant admixture of the higher lying 
states. However, the gates must be switched smoothly 
enough for the evolution to be adiabatic with respect to 
the higher lying states, as will be seen later from the 
numerical analysis. 

The energy spectrum in Fig.|3|3 is periodic in the mag- 
netic field in the usual AB sense. Since tp relates to one 
quarter of the phase on the entire outer loop, this means 
that with every Aip — 2-7r/4 one additional flux quantum 
enters or leaves the cross-sectional area of the array. This 
is seen for example in the splitting of the singlet which 
opens and closes with a period Atp = 2tt/A. The exact 
period of the system, however, is Atp = it. Note, that 
this is not because of the usual t — > — t symmetry which 
does not hold here because of the diagonal cross link in 
the array shown in Fig. Instead, it can be related 
to changing the sign in both basis states of the qu2LS. 
This is easily seen for tp = it from Fig.[5jD by considering 
0}^ h = 1 13) = — 1 31) and dropping all the arrows shown. 

Now the essential effect of the magnetic field is that it 
allows one to close the gap in the singlet qu2LS while at 
the same time it opens a gap in the triplet qu2LS (at fixed 
S z ). The smallest magnetic field where this happens is 
at tp = ipo — 0.286 7r, indicated by the arrow in Fig. [2Jd. 
The important consequence of tuning the magnetic field 
to tp — ► tpo is that the charge qubit can be held frozen in 
its state by also having V g — (see below). 

Figure 13; and d show the singlet ground state configu- 
ration and its two spatially distinct basis states, respec- 
tively. Note that in order to have (close to) degenerate 
groundstate qu2LS, there must exist a basis representa- 
tion that is spatially complementary 9 , in agreement with 
what is shown in panel (d). 

In the context of quantum computation, any linear su- 
perposition of the two states forming the qubit must be 
possible^ The problem is conveniently mapped into a 
pseudo spin Hamiltonian with its equivalent 2-level sys- 
tem. The general qubit state \tp) = (ci,C2) when written 
as a density matrix p — \tp) (tp\ = | (1 + r^a) defines the 
3D Bloch vector rj, which then is used as the represen- 
tation of the qubit stated Within this frame, an arbi- 
trary qubit operation translates into a rotation of the 
Bloch vector (pseudo spin) anywhere in its 3D domain. 
Two distinct rotations are sufficient to do so. The rota- 
tions in pseudo spin language are then generated by the 
Pauli matrices o~{ x ,y,z} and the requirement of two dis- 
tinct rotations translates into two distinct quantum gates 
(qugates) that can be built from the Pauli matrices. 

For the charge qubit encoded in the 2x2 array with 
the basis states as shown in Fig. the physical quan- 
tum gates are now as follows: the asymmetrically applied 
gate voltage (see Fig.^i) only drags apart the potentials 
of the two qubit basis states in Fig. and thus rep- 
resents the c z -gate, also referred to as V-g&te. On the 
other hand, the singlet gap or anticrossing can be re- 
lated to a real off-diagonal element in the 2D pseudospin 
Hamiltonian which can be tuned down to exactly zero 
by a magnetic field as shown above. Thus this is referred 
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FIG. 3: Energy spectra for the qu2LS of the 2x2 qudot sys- 
tem together with a few higher lying states for singlet and 
triplet states (\S, S z ) = |0, 0) and \S, S z ) = \l,m= -1,0,1} 
respectively), (a) Energy spectrum vs. asymmetrically ap- 
plied gate voltage, Vg = V gi = — V g2 . The doubly occupied 
states lie about lmeV higher in energy (outside figure) and 
therefore have negligible influence. The inset shows a closeup 
of the (anti)crossing in the qu2LS. (b) Energy spectrum vs. 
uniform external magnetic field perpendicular to the array ex- 
pressed through the phase in t = \t\ e tlp . The initial singlet an- 
ticrossing at ip = is completely closed for if = Lpo — 0.286 n, 
indicated by the arrow in panel (b), while at the same time 
the triplet levels show a pronounced anticrossing. (c) Singlet 
ground state probability distribution over the 2x2 array. This 
state is a symmetric combination of the basis states shown in 
panel (d): Probability distribution of the basis states of the 
singlet qu2LS labelled \0} qb 
to find spin up or spin down 



and \l) qb with equal probability 

IvWI 2 = IVul 2 - 



to as the cr^-gate or i?-gate. Together, the two physical 
qugates introduced can be utilized to generate arbitrary 
rotations of the Bloch vector and thus to construct ar- 
bitrary qugates for the qu2LS. Also, since both of the 
qugates can be turned off completely by setting V g = 
and turning on a specified magnetic field, this allows to 
freeze the qu2LS in any arbitrary state at any time. 



A. Splitting due to exchange energy 

We want to analyze now the reasons for the closing 
and opening of the gap in the qu2LS, as a way to provide 



us with insights into the nature of these states. The 
Hamiltonian in Eq. is written out explicitly as a six- 
dimensional matrix in Eq. © for two electrons on the 2 x 
2 array where, for simplicity, states of double occupancy 
are neglected. 
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The sign in — (+)t refers to the singlet (triplet) spatial 
Hamiltonian, respectively, with the complex tunneling t 
as in Eq. (0). The top row and the left-most column of 
Eq. indicate the two-particle basis states chosen as 
shown in Fig. \Bp. With t — 0, the qu2LS grouped in the 
upper left 2x2 block of the H matrix is degenerate and 
all intermediate states are split off by the same energy 
Ao = £2 — £i due to symmetry. Here, £i and £2 are the 
overall diagonal contributions arising from the Coulomb 
interaction at V g = 0. 

The Hamiltonian in Eq. © can be diagonalized an- 
alytically. However, a perturbative approach provides 
explicit insights on the reasons for the splitting due to 
exchange of otherwise degenerate states. Furthermore, 
the setup proves sufficiently simple to allow the complex 
sum to all orders over all possible histories in Hilbcrt 
space within the Feshbach formalism. The result consis- 
tently agrees with the analytical solution to the problem. 

In order to obtain an estimate for the splitting in the 
ground state set, the Feshbach formalism provides an ef- 
fective Hamiltonian for the qu2LS (V space) coupled to 
higher lying states (Q space, see App. 0). For conve- 
nience, the qu2LS basis is written as V = {|13} , |24)} = 

{\Q) qh ,\l) qh \ = {0, l} gb . The remaining intermediate 

states in Fig. |2b (Q space) form a ring topology where 

Note 
but 



every node is linked to the |0) 9b and |l) o6 states 



Iqb 



that there is no direct transition from \Q) qb to |1) 



Iqb 



one has to proceed through at least one of the interme- 
diate states with an energy cost of A = £2 — £i- 

The effective 2D Hamiltonian for the qu2LS in the ab- 
sence of a magnetic field is now constructed as follows: 
the matrix element (H e ff)^ with i,j = {0,l} gb is the 
sum over all possible paths that start in state i, immedi- 
ately proceed to intermediate states Q, and only in the 
final step come back to state j. The number of possible 
paths constructed in this manner for a total of n steps is 
S n = 2™ with n > 2 since there are 4 possibilities 
to go from |i) to the intermediate states, then 2 possibil- 
ities of which way to go in the ring at each of the n — 2 
intermediate steps, and 1 choice left to finally leave the 
ring and go to state \j). Furthermore, for n > 2, exactly 
half of these paths include particle exchange, i.e. have 
an odd number of dashed lines in Fig. ^jp, and thus for 
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the case of the triplet states cancel each other to zero. 
The underlying reason for this is the spatial C^ v symme- 
try of the 2x2 setup which has two mirror symmetry 
planes perpendicular to the array. Therefore, for every 
path starting from {0, l} qb and ending in {1, 0} qb has a 
mirrored counterpart where the particles are exchanged 
in the final state as compared to the first path. 

The case n = 2 needs separate consideration. For the 
diagonal elements in H e g the same step is taken twice, 
back and forth, and thus the relative sign in t does not 
matter. Therefore the triplet states have an n = 2 contri- 
bution in the diagonal. Putting all these pieces together, 
yields the matrix elements for the effective Hamiltonian 
which, for example, in case of the singlet states are 



(i\ H eff \j) S = £ l S ij 



(-2tf 



oj - e 2 + 2t 



(~2ty ( -2t 

e 1 6 ij + Eg M , 



(4) 



and thus 



H 



singlet 
eff 



H 



triplet 




where for the triplet state only the paths with n = 2 
contribute. The eigenstates for the qu2LS are now ob- 
tained from the nonlinear eigensystem H e g (lj) = 
uj | -0) where is restricted to the 2D ground space. 
jjtnp et j g g £jjj diagonal, and therefore the triplet states 
do not mix with each other, but are just shifted lower by 

-o m 



^ - V^) + ^ = + u \T) > Wltn A ° - e 2 - £l . 

The singlet states, however, rearrange to symmetric and 
antisymmetric combinations of the original basis. One 
of the eigenstates remains at the original eigenenergy 
uj = ei, while the second one is lowered by S given as 



(6) 

and thus forms the ground state of the system for finite 
t. The analytical solutions are consistent with the nu- 
merical diagonalization of the Hamiltonian in Eq. 
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FIG. 4: Energy splitting in the qu2LS for the system of 
Eq. |3] in dependence of the magnetic field for singlet (red) 
and triplet states (blue). The splitting is shown in units of 
Ao = £2 — si, namely the separation of the qu2LS from the 
remaining Hilbert space. The dashed black line is the result 
of the lowest order Feshbach analysis, Eq. J7J. 



since they correspond to transitions through the diago- 
nal of the array in Fig. |2t (there are no arrows on the 
transitions between intermediate states in Fig. [2b) . 

The paths can be summed up similarly. For simplicity, 
however, we give only the contribution to lowest order in 
t for the effective Hamiltonian. The denominator uj — e 2 
is then replaced by E\ — e 2 = — Ao, and so the lowest 
order contribution to the self energy term is 



^singlet s _ / 1 COS 2ip 

Aq I cos 2ip 1 , 



J eff 



+ 0(t 3 ) (7a) 



^triplet 
J zff 



Ao \-isin2y? 1 



Comparing this with the last terms in Eqs. (JSJ, the ef- 
fect of the external magnetic field is obvious. With in- 
creasing ip the singlet splitting can be reduced down to 
zero, while simultaneously a comparable gap opens in the 
triplet states, in agreement with the numerical data for 
the full 2x2 system in Fig. |3t>- The effect of the mag- 
netic field on the singlet (triplet) states is that of a a x 
(cry) gate, and thus clearly provides the necessary second 
quantum gate for single qubit operation. Figure ^com- 
pares the exact numerical results of the Hamiltonian in 
Eq. © with above lowest order perturbative approach for 
the splitting in the qu2LS. The lowest order contribution 
in the Feshbach formalism already provides an excellent 
approximation. 



B. Effect of a magnetic field 

We apply again the Feshbach formalism to the system 
in Fig. [2b, but now including the complex phases indi- 
cated by the arrows in that figure. Note that the phases 
due to the magnetic field affect only the first and last 
step in each path, while intermediate state transitions 
remain unaltered by the presence of the magnetic field 



C. Numerical qubit dynamic 

The time evolution of the 2x2 qudot array is studied 
numerically for the singlet state under the action of the 
V-gate (a z ) and the .B-gate (a x ). Note that here the 
.B-gate is considered non-active when a magnetic field 
is tuned to the phase tp = tpo indicated by the arrow 
in Fig. |2h>, while the gate is considered active when the 
magnetic field is turned off (ip = 0). In this sense, with 
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none of the two gates applied, the system is static since 
the singlet states are degenerate. Typical time dynamics 
data is shown in Fig. Starting in the singlet ground 
state (Bloch vector f% — +x), the system is static at 
t = 0. The following U-gate rotates this state around 
the z axis by 450°, ending in fb = +y, where the system 
is stalled for a short time interval until the B-gate is 
activated at t = 1.6ns. Now the S-gate rotates r% around 
the i-axis again by 450°, leaving the state in the +z 
state where the system is stalled again at t = 3ns, now 
being in the central region on the time axis in figure 
Fig. [SJd. Since ±z corresponds to the basis states of the 
qubit, the charge distribution equals the \l) qb state in 
Fig. OJi as can be seen by the snapshots shown along the 
time evolution in between panels (a) and (b) in Fig. 
After another B-gate of the same duration, the system is 
rotated another 450° around the i-axis leaving the qubit 
in the — y state at t = 4.8ns. When finally a (— U)-gate 
is applied, the qubit is rotated by —450° around the z- 
axis and the system is left in the —x configuration at 
t = 6.2ns. The time evolution of the Bloch vector in 
this whole process sweeps two grand circles in the Bloch 
sphere, as shown in the inset of panel (b) . 

The numerical data shown in Fig. confirms the pre- 
vious analysis. The main point is that the qubit can be 
placed into any state by applying appropriate magnetic 
field and gate voltages to the system, yielding full control 
of the qubit as desired. 



IV. CONCLUSIONS 

Singlet and triplet states have been studied on a planar 
array of quantum dots. With respect to qubits, a two-fold 
nearly degenerate ground state pair was constructed and 
the splitting of this two-level system was explained and 
estimated using the Feshbach formalism. Furthermore 
these states were shown to represent full single qubit op- 
erations encoded in the charge states. The AB flux given 
by a uniform magnetic field provides the required second 
quantum gate by generating a dynamical phase in the 
wave function of the two electron qubit state. This can 
be used to effectively turn off the qubit dynamics at will. 
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FIG. 5: Time evolution and control of the singlet qu2LS 
on the 2x2 array, (a) Time evolution of the state occu- 
pancy with respect to the qubit basis \0) qb and \1) b (see 

Fig-Eli)- (b) Time evolution of the site occupancy |(c^Ci)| 2 = 

\(ip \cfci\ tp)\ 2 ■ The square panels in between panels (a) and 
(b) show snapshots of the charge distribution on the array 
at the times indicated either towards panel (a) or panel (b). 
The inset in panel (b) shows the time evolution of the qubit 
in Bloch sphere representation, (c) Time dependence of the 
voltage gates (black) and the magnetic field expressed through 
Re (t) and Im (t) (red lines) where Abs(t) is kept constant. 
The time constant for rise and fall time of the gate voltages 
was chosen as Ty = 0.658 ps while for the tunneling the con- 
siderably longer t v = 100 • Ty = 65.8 ps was used out of 
adiabatic purposes with respect to the higher lying states. 



infinite) system to the Hamiltonian of a small (finite) sub- 
system which should be energetically well separated from 
the remainder of the spaced Given a finite subspace V of 
the total Hilbert space Tt with its complement Q, such 
that V + Q = H, the projections into these spaces are 

P = £^I*X*I and Q = J2k ir \k)(k\ = 1 - P, re- 
spectively. The stationary Schrodinger equation when 
projected into the V and Q spaces becomes 



Hpp Hpq 
H QP Hqq 




= E 




(Al) 



APPENDIX A: FESHBACH FORMALISM 

The Feshbach formalism provides an efficient proce- 
dure to reduce the full Hamiltonian of a large (possibly 



with the projections defined as Hpq = PHQ, \ipp) = 
P \ip), and similarly for the remaining ones. By assump- 
tion, V has finite dimension, thus Hpp is also finite. \iPq) 
can be formally eliminated from Eq. IjAlf) and the result 



7 



is H^ff \ipp) = E \tpp) with 

H eff = HpP + Hp Q F _ tt H QP> ( A2 ) 



=Sc3Q(S) 

with H^ff describing the effective Hamiltonian in space 
V . H^ff is thus the sum of the unperturbed matrix ele- 
ments Hpp and the self-energy contribution Sqq. Note 
that Eq. HA2fl is still exact and thus must be nonlinear in 
E in order to represent the entire eigenspectrum of the 
system. The non-linearity is manifested in the E depen- 
dence of Sqq. It ensures a good approximation while it 
does not result in further complications here, since for 
the two-dimensional space V discussed in this paper, the 
resulting equations can be conveniently solved analyti- 
cally. 

If the full Hamiltonian H did not couple the V and 
Q spaces, then with Hpq = the second term Sqq in 
Eq. l|A2(l vanishes consistently. Note also that eliminat- 
ing \i[>q) from Eq. I|A1|) eliminates the coefficients of \i/)q) 



so that this procedure systematically eliminates variables 
on the large scale. Eq. I|A2[1 is a formal solution since the 
self-energy Sqq (E) is not known. But there are different 
ways to approximate £qq using perturbative expansions 
naturally suggested by its definition. With the identity 



cj - (H + V) QQ uj - H 0:QQ ^ V w - h o,qq J 

(A3) 

a very useful generalization of the Brillouin-Wigner for- 
malism to more than one, but still a finite number of 
states, is obtained^ This approach is essentially a path 
formulation in Hilbert space in the sense that all possible 
path histories through Hilbert space are taken, starting in 
V, proceeding directly into Q (Hpq in Ea. lA2|) and only 
in the final step returning back to V (Hqp in Ea. IA2J1 . 
In addition, every intermediate step is weighted by the 
propagator terms (k\ (E — Hqq) \k) = (E — £&) 1 
fEa. lA3|) . The lowest order contributions are then equiva- 
lent to the shortest histories through Q space with lowest 
cost in energy. 
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